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COST  BAS  CD  ACCEPTANCE  S.WLINC  PLANS 
AND  PROCESS  CONTROL  SCHEMES 


Douglas  C.  Montgomery 

School  of  Industrial  and  Systems  Engineering 
Georgia  Institute  of  Technology 
Atlanta,  Georgia  30332 


ABSTRACT 

A  review  of  procedures  for  designing  acceptance 
sampling  plans  and  process  control  schemes  according 
to  economic  criteria  Is  given.  The  use  of  nonstt- 
tlsclcal  criteria  to  design  these  procedures  has 
bcea  the  subject  of  considerable  research,  and  a 
number  of  different  process  models  and  procedures 
have  been  developed.  A  comparison  of  the  major  ap¬ 
proaches  to  model  formulation  and  design,  and  a  dis¬ 
cussion  of  ths  praccicsl  implementation  of  these 
techniques  is  the  primary  focus  of  this  paper. 

INTRODUCTION 

Acceptance  sampling  and  process  control  are  the 
primary  statistical  techniques  used  in  a  quality 
assurahca  program.  These  techniques  have  found  wide 
application  in  industry,  particularly  in  manufactur¬ 
ing,  in  the  following  arena:  Incoming  material 
laspactlon,  surveillance  of  the  production  proceas, 
estimation  of  lot  or  procesa  characteristics,  process 
capability  analysis,  and  finished  product  quality 
auditing. 

Traditionally,  acceptance  sampling  and  process 
Control  schemes  are  designed  with  respect  to  etetis- 
tlcal  criteria.  For  example,  one  may  choose  the  sam¬ 
ple  sirs  and  acceptance  number  for  a  single-sampling 
plan  far  attributes  so  that  tha  operating  character¬ 
istic  curve  passea  through  (or  near)  two  points 
specifically  selected  to  give  certain  probabilities 
of  lot  rejection  et  specified  levels  of  lot  or  pro¬ 
cess  quality.  Similarly,  one  say  design  s  control 
chart  so  that  the  power  of  the  chart  to  detect  a 
particular  shift  in  process  quality  and  the  probabi¬ 
lity  of  false  alarms  are  equal  to  specified  values. 
While  tha  traditional  approach  often  produces  accsp- 
table  results,  it  la  also  possible  to  design  accep¬ 
tance  sampling  and  process  control  schemas  with  re¬ 
spect  to  economic  criteria.  This  has  considerable 
intuitive  appeal  since  there  has  been  an  increasing 
emphasis  on  quality  costs  In  recent  years.  Further¬ 
more,  tha  use  of  these  techniques  have  direct  econo- 
mlc  consequences  in  that  one  la  balancing  the  costs 
associated  with  sampling,  tasting  and  process  sur¬ 
veillance  against  internal  and  external  failure 
costs.  Since  the  design  of  the  procedure  affects 
these  costa,  it  la  logical  to  consider  this  design 
from  an  economic  viewpoint. 


During  the  last  ZO  years,  considerable  re¬ 
search  has  been  devoted  to  the  use  of  economic  cri¬ 
teria  to  design  acceptance  sampling  plana  and  pro¬ 
ceas  control  schemes.  This  paper  gives  an  overview 
of  the  major  developments  and  approaches.  Some 
comparison  of  the  different  approaches  to  model 
formulation  and  a  discussion  of  the  practical  imple¬ 
mentation  of  these  techniques  is  also  given. 

ECONOMIC  MODELING  AND  PROCESS  CONTROL 

The  Shevhart  control  chart  is  probably  the  most 
widely  used  process  surveillance  device.  To  design 
a  control  chart,  we  oust  choose  the  sample  size  (n), 
the  control  limit  (k,  the  multiple  of  o),  sod  the 
interval  between  samples  (h  hours) .  Despite  tha 
non-optimality  of  fixed  sample  size,  fixed  sampling 
Interval  procedures,  the  Shevhart  control  chart  has 
gained  widespread  use  because  of  its  flexibility, 
simplicity  of  administration  and  tha  additional  in¬ 
formation  about  procesa  performance  often  contained 
in  the  pattern  of  points  plotted  on  the  chart. 

Contiderabla  effort  haa  been  devoted  to  deve¬ 
loping  economic  models  of  Shevhart  control  charts. 
These  models  ususlly  assume  Chat  the  process  is 
characterized  by  one  in-control  state  Chat  represdnta 
the  mean  of  the  quality  characteristic  whan  no 
assignable  causes  are  present,  and  s  2  1  out-of- 
control  states.  The  probability  modal  that  governs 
the  transitions  between  these  s  +  1  statee  is  called 
the  process  failure  mechanism. 

Costs  and  Measures  of  Ef fsctlvenass 

Three  cstegorles  of  costs  are  usually  consid¬ 
ered  in  the  dsvelopswnt  of  economic  models  for 
Sbewhart  control  charts:  tba  cost  of  sampling  and 
tasting,  usually  of  the  form  aj  +  ajn,  thn  costa  of 
investigating  and  possibly  correcting  action  signals, 
•end  tha  coat  of  producing  defective  ltams  (internal 
and  axtsrnsl  failure  costs) .  These  costa  are  than 
comblasd  to  form  a  total  cost  par  unit  time  function. 
Ths  general  approach  is  to  define  e  cycle  ee  Che 
length  ef  time  T  during  which  the  proceee  begins 
operating  In  the  ln-control  etete  end  eventually  re¬ 
turns  to  ths  ln-control  state  following  e  proceee 
edjuataenc. 
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If  C  la  the  cost  Incurred  during  a  cycle,  then 
Che  expected  coat  per  unlc  time  la 

EU)  -  E(C)/E(T)  (1) 

Since  E(C)  and  E(T)  are  functions  of  Che  design  par* 
enters  of  the  control  chart,  numerical  optimization 
(direct  search)  techniques  can  be  employed  to  find 
the  parameter  values  that  minimize  the  expected  cost 
per  unit  time. 

Basic  Process  Models 

The  fundenental  process  model  often  used  in 
ecoooalc  design  of  the  Shewhart  control  c..art  Is  due 
to  IXrncan  (7J.  His  wort  was  Che  first  ca  deal  with 
a  fully-economic  approach  to  control  chart  design 
and  to  explicitly  consider  the  optimization  problem. 
The  Ikincan  model  assumes  a  single  assignable  cause 
(thus  a-1)  and  Chat  transitions  batvean  the  in- 
control  and  out-of-control  states  occur  according 
to  a  Poisson  process  with  intensity  X  occurrencss 
per  unit  time.  Thus  the  length  of  time  the  process 
remains  in  control,  given  that  it  starts  In  that 
state.  Is  an  exponential  random  variable  with  mean 
1/1.  This  is  equivalent  to  a  process  failure  mech¬ 
anism  with  a  uniform  hazard  function. 

The  cycle  conelsca  of  four  perlode:  (1)  the 
in -control  period,  (2)  the  out-of-control  period,  (3) 
the  time  to  take  a  sample  and  lnterprat  the  resulte 
and  (4)  the  time  to  find  the  aaslgnable  cause.  As¬ 
suming  that  the  process  continues  to  run  during 
searches,  the  expected  cycle  length  le 

E(T)  -  1/1  +  h/U-B)  -  T  +  gn  +  D  (2) 

where  h  la  the  sampling  interval  (in  hours) ,1-8 
is  the  power  of  the  chart  for  detecting  a  specified 
out-of-control  state,  n  Is  the  sample  site,  g  is 
the  tine  required  to  take  and  Interpret  e  sample,  D 
la  the  expected  time  to  find  an  assignable  cause  and 

1  (J+l)h  (J+l)h 

T  •J  l(t-jh)dt  /  f  «*“  Idt 

1*  jh 

|  -  (1  -  <1+Ah)e“*h]  /  l(l-e'Xh)  (3) 

la  the  expected  time  of  occurrence  of  the  shift, 
glvm  Chat  It  occurs  betweeen  the  and  (J+l)— 

sample*. 

Let  Vo  and  V],  be  the  net  income  per  hour  of 
operation  in  the  in-control  and  out-of-control 
'State*,  respectively .  The  cost*  of  investigating 
real  and  false  alarm*  are  *3  and  si,  respectively, 
lines. there  are  E(I)/h  samples  per  cycle  and  a*-"1/ 
(l-*~  false  alarm*  par  cycle,  on  th*  average,  the 
expected  net  income  par  cycle  is 

1(C)  •  V0(l/X)  +  Ux*m2n>  «(T)/h 
♦  Vl[h/(1-B)  -  x  ♦  gn  +  0] 

-  *j  -  *j«*"*h  /  (l-e*ih)  (4) 


Dividing  (4)  by  (2)  gives  the  expected  net  in¬ 
come  per  hut.r,  say 

E(I)  -  V0  -  E(L)  (5) 

where 

E(L)  *  (»j+a2n)/h 

+  Je^[h/(1-B)  -  t  +  gn  ►  D] 

.  ,  *  -Xh  ,  ,,  -Xh,( 

♦  *3  +  Sjoe  /  (1-e  ){ 

♦  [1/X  ♦  h/(l-B)  -  t  +  gn  ♦  D]  (6) 

and  a^  ■  Vq  -  Vy  ia  the  penalty  cost  (per  hour)  re¬ 
sulting  from  production  in  the  out-of  control  state. 
Maximizing  Ed)  In  (5)  is  equivalent  to  minimizing 
E(L)  in  (6).  There  are  three  decision  parameters, 
n,  k  sad  h  (the  function  E(L)  depends  on  the  width 
of  the  control  limits  k  through  a  end  6) ;  end  nine 
user-supplied  constants  that  describe  the  production 
process,  4^,  42,  *3.  *3.  «4,X,  g,  D,  and  6  tha  mag¬ 
nitude  of  the  process  shift  (which  la  Involved  in 
determining  o  and  B)> 

While  this  bsalc  process  model  could  be  used 
in  many  situations,  there  are  a  number  of  variations 
that  are  potentially  important.  The  first  of  these 
is  the  incorporation  of  more  than  one  out-of-control 
state.  Duncan  [8]  haa  generalized  the  cost  model 
above  to  a  situation  where  there  ere  s  out-of -control 
states  (aaslgnable  causes).  Knappanberger  and  Cran- 
dage  [13]  also  present  a  model  capable  of  treating 
several  out-of-control  state*.  The  Duncan  model 
appears  somewhat  more  realistic  in  that  it  employ* 
different  search  costs  for  different  aaslgnable 
causes,  while  the  Knappanberger  and  Crandege  model 
uses  an  average  or  composite  coat.  On  the  other 
hand,  thq  Kneppenbcrgsr  and  Grandag*  modal  allow* 
continued  process  deterioration  beyond  the  initial 
shift,  while  th*  Duncan  modal  allow*  only  on*  or  two 
shifts  to  occur.  Both  models  ars  considerably  more 
sophisticated  in  terms  of  their  data  requirements 
than  (6)  above,  furthermore,  there  is  considerable 
evidence  Chat  a  single  assignable  cause  model  that 
matches  the  true  multi-state  syatam  in  important 
weye  will  be  a  satisfactory  approximation.  Specifi¬ 
cally,  if  we  let  the  single  out-of-control  state  be 
a  weighted  eysrage  of  th*  s  out-of-control  states 
with  th*  weights  chosen  proportionally  to  tha  proba¬ 
bility  of  occurrence  of  th*  state,  excellent  result* 
at*  obtained.  In  prsctic*,  there  is  probably  lltcle 
reason  eo  use  a  multiple  out-of-control  state  model 
in  most  situation*. 

Tha  assumption  of  a  constant  hazard  function 
for  th*  proc***  failure  mechanism  (an  exponential 
distribution  for  tha  time  between  shifts)  le  criti¬ 
cal.  If  tha  occurrence  of  assignable  causes  can  be 
thought  of  as  random  "shock*"  disturbing  the  systmt, 
then  the  probability  of  a  thift  occurring  within  any 
small  interval  of  time  ia  proportional  to  th*  length 
of  th*  interval,  and  th*  expos  at  i*l  model  is  prob¬ 
ably  appropriate.  However,  if  assignable  causes 
occur  as  a  result  of  tha  cumulative  effects  of  vibra¬ 
tion  or  hast,  improper  set-up,  or  excooelvs  otroooes 
during  start -up,  a*  Is  th*  coso  with  nany  production 


processes,  then  the  exponential  model  may  be  entirely 
Inappropriate.  Distributions  with  either  increasing 
or  decreasing  hazard  functions  may  be  more  reasonable. 
Misspecif  icatlon  of  this  aspect'  of  the  process  model 
may  have  very  serious  economic  Impact.  Unfortunately, 
It  Is  not  straightforward  to  Incorporate  other  process 
failure  mechanisms  in  (6).  A  discrete-time,  single¬ 
cause  model  that  allows  free  choice  of  the  process 
failure  mechanism  has  been  developed  by  Baker  [2], 
but  as  it  assumes  that  a  sample  Is  taken  at  the  end 
of  each  period,  the  optimum  Intersample  interval  can¬ 
not  be  explicitly  determined. 

The  assumptions  In  (6)  that  the  process  contin¬ 
ues  in  operation  during  the  search  for  ar  assignable 
cause  end  that  the  co-t  of  repair  la  not  Included  are 
unrealistic  for  many  processes.  Fortunately,  (6)  may 
be  easily  modified  to  consider  these  possibilities. 

If  Dq  and  D}  are  the  expected  search  times  for  false 
and  real  alarms,  respectively,  and  If  A  is  the  cost 
of  repair,  then  assuming  that  the  process  la  shut 
down  during  the  search  leads  to  an  expected  net  In¬ 
come  per  unit  time  of 

E(I)  -  V0  -  EO.)  (7) 

where 

E(L)  -  J(a1+a2n)U/X  +  h/(l-e)  -  t]/h 

+  «3  ♦  A  +  [«m~Ah/(l-e_Xh)1(V0D0  +  aj) 

+  V0Dq  +  a4Ih/U-«  -  Tl| 

♦  [1/x  ♦  h/U-B)  -  T  +  ca>0a“Xh/ 

(l-*‘xh)  +  Dx]  <g) 

Tor  a  detailed  development  of  this  model,  see  Mont¬ 
gomery  [U],  Rota  that  given  alt  a2,  *3,  aj,  A,  •*, 

X,  Vo*  ®0»  ®i*  the  magnitude  of  the  process  shift, 
(8)  may  6a  easily  optimised  for  n,  h,  and  h  using 
direct  search  net hod a. 

Malar  Applications 

Suppose  that  the  quality  characteristic  of  In¬ 
terest  la  represented  by  the  parameter  8,  mad  that  a 
represents  the  eample  statistic  corresponding  to  9. 
There  ere  two  possible  values  for  9;  6  -  60  corres¬ 
ponding  to  the  la-ceotrol  state  and  8  •  82  which 
represents  the  out -of -control  state.  Then  the  pro¬ 
bability  of  a  false  alarm  is  : 

a  -  1  -  TOCL  i  x  i  UCL  |  8  -  BQ)  CD 

and  the  power  of  the  control  chart  it 

1  -  8  -  1  -  TOO.  «  *  S  VCL  |  •  -  9t)  (10) 

whore  ICl  and  UCL  sad  thp  lower  and  upper  control 
limits  am  the  Shewhert  control  chart.  This  by  spe¬ 
cifying  the  quality  characteristic ,  the  sample  sta¬ 
tistic,  sod  the  relay ant  sampling  distribution,  ths 
coot  models  of  tho  previous  saction  could  ba  applied 


to  any  type  of  Shewhart  control  chart.  For  example, 
if  the  quality  characteristic  of  Interest  Is  a  mea¬ 
surement  described  by  the  process  mean,  then  the 
x-ct.art  would  be  used;  thus  equations  (9)  and  (10) 
become 

a  *  28(-k) 
and 

1  -  8  -  *(«v£-k)  +  »<-«v£-k) 

respectively,  where  *(t)  denotes  the  standard  normal 
distribution  function. 

The  major  applications  of  cost  models  suchas 
those  In  the  previous  section  have  been  to  the  x- 
chart,  the  x/R  chart  combination,  and  the  fraction 
defectlva  or  p-chart.  Some  work  has  also  focused  on 
procedures  for  the  simultaneous  control  of  several 
related  quality  characteristics,  and  the  treatment 
of  non-Shewhart  control  charts  such  as  the  cumulative 
sum  control  chart  and  the  use  of  warning  limits  on 
tha  x- chart .  Almost  all  reported  studies  assume  that 
the  procaaa  failure  mechanism  la  exponential.  As 
noted  earlier,  this  la  a  critical  and  sometimes  un¬ 
warranted  assumption.  For  an  extensive  review  of  the 
literature  in  this  area,  see  Montgomery  (14]. 

Most  of  ths  papers  cited  In  114]  contain  numer¬ 
ical  studies  and  usually,  a  sensitivity  analysis. 

From  theca  results,  it  la  possible  to  draw  certain 
general  conclusion*  about  control  chart  design: 

1.  Tha  optimum  sample  slsa  la  largely  determined  by 
ths  size  of  tha  shift,  with  smaller  shifts  re¬ 
quiring  larger  samples.  Small  shifts  may  require 
very  large  samples,  possibly  as  large  as  a-40  or 
mors  on  an  x  chart. 

2.  Changes  in  a2  and  a2,  tbs  fixed  and  variable  costs 
of  sampling,  affect  all  three  design  parameters. 
Increasing  ths  fixed  cost  a.  increases  tha  Inter¬ 
val  Between  samples  and  leads  usually  to  slightly 
larger  samples,  while  Increasing  tha  variable  cost 
aj  usually  results  la  small,  infrequent  samples, 
hit  narrow  control  limits. 

3.  The  penalty  cost  for  out -of -control  production 
a*  safely  Influence*  tha  interval  between  asm- 
plea  b.  _  Larger  value*  of  a*  imply  emaller  opti¬ 
mum  values  of  h..  A  similar  affect  is  observed 
by  increasing  X,  the  mean  number  of  process 
shifts  per  hour. 

■A.  Tha  search  coats  aj  and  a«  mainly  affect  the  con¬ 
trol  limits,  end  have  s  slight  affect  00  sample 
else.  Larger  values  of  03  sad  aj  result  in  wi¬ 
der  control  Its  It  a  and  larger  sample  alee  a,  as 
increasing  ssareh  tests  imply  that  fewer  false 
slaves  ar*  desirable. 

3.  Ths  optimum  economic  control  chart  design  Is 
relatively  Insensitive  to  estimates  of  tho  coot 
parameters.  Ths  cost  surfaces  are  usually  flat 
near  the  optimum,  and  are  usually  stsopar  soar 
tha  origin ,  so  that  It  is  hotter  to  overestimate 
the  control  chart  donlgn  parameters  than  te  un¬ 
derestimate  them.  Ths  optimum  design  is  moot 


« 


sensitive  to  errors  in  ostimating  the  In- 
control  and  ouc-of-oontrol  states  (or  the  Mag¬ 
nitude  of  the  process  shift) . 


Once  the  opting!  n  and  k  are  determined,  the  inter¬ 
val  between  samples  could  be  obtained  from 


this  last  finding  has  important  practical  impli¬ 
cations,  for  it  is  often  difficult  to  precisely  esti¬ 
mate  costs.  However,  the  mean  time  between  shifts 
and  the  magnitude  of  the  process  shifts  are  usually 
mota  easily  determined  from  process  performance 
data  or  from  the  engineer's  knowledge  of  the  opera¬ 
ting  environment.  Note  that  one  should  exercise 
caution  in  using  arbitrarily  designed  control  charts 
and  empirical  "rulaa  of  thumb",  such  as  n»5,  k*3  and 
1»*1  for  the  x-chart.  In  some  cases,  particularly 
those  with  small  shifts  and  large  penalty  costs 
for  production  in  the  out-of-control  state,  signifi¬ 
cant  economic  penalties  may  result. 

A  Seml-Economlc  Approach 

Instead  of  the  fully-economic  modeling  approach 
cxitllnad  above,  it  may  be  desirable  in  some  cases  to 
use  a  semi-economic  scheme  that  blends  both  economic 
and  statistical  criteria  directly.  This  could  be 
useful  in  situations  where  the  analyst  la  unwilling 
or  unable  to  estimate  all  of  the  fully-economic 
model's  coat  parameters,  or  where  he  wishes  a  com¬ 
promise  between  the  fully-economic  and  purely  etatia- 
tical  approaches. 

Oat  such  semi- economic  scheme  would  be  to  find 
the  values  of  n,  k,  and  h  that  minimize  the  expected 
sampling  costs  per  cycle,  but  that  also  have  desir¬ 
able  average  run  lengths  for  the  control  chart  in 


both  the  in-control  and  aut-of -control  states. 

Can- 

sequeatly, 

we  would  like  to  choose  n,  k,  and  h 

SO  AS 

to 

minimize  Z  *  (aj+«2n)  E(T)/b 

Ul) 

- 

maximize  AttQ  «  1/a 

(12) 

minimise  Att^  ■  1/U.-8) 

(13) 

h  • 


i  \ 

aa3  +  ai  +  a2o  > 
la^Kl-e)-1  -  0.5}  I 


(15) 


This  is  an  approximation  for  the  optimum  sampling 
interval  h  given  n  and  k  suggested  by  Duncan  [7] 
and  Chiu  and  Wetherlll  [4], 

These  formulations  of  the  optimal  control  chsrt 
design  problem  have  not  been  extensively  investi¬ 
gated.  It  would  be  of  interest  to  discover  how  they 
compare  with  the  fully-economic  solutions,  and  with 
the  standard  designs  often  suggested  In  the  litera¬ 
ture. 


ECONOMIC  ASPECTS  OF  ACCEPTANCE  SAMPLXNC 

There  are  various  types  of  acceptance  sampling 
schemes;  single  sampling,  double  sampling,  sequen¬ 
tial  sampling,  continuous  sampling,  and  ao  forth. 

The  quality  characteristic  being  Inspected  may  be 
either  an  attribute  or  a  variable.  For  various 
reasons,  more  attention  has  been  given  to  attrlbutaa 
sampling ,  and,  in  particular,  single  sampling  plans 
for  attributes  have  been  studied  extensively.  He 
now  give  some  of  the  key  results  pertaining  to  the 
economic  design  of  these  plana  and  briefly  survey 
some  of  the  other  related  work. 

An  attribute  single  sampling  plan  is  indexed  by 
three  numbers:  the  lot  size  N,  the  sample  site  n, 
end  the  acceptance  number  c.  If,  In  a  random  sam¬ 
ple  of  size  n,  c  or  fever  defectives  ars  found,  tba 
lot  is  accepted,  while  if  more  than  c  defectives  ate 
found  the  lot  is  rejected.  The  probability  of  lot 
acceptance  Is 

c  /  K0\  /H-N8  \  //  N  \ 

?<*>  -  ^  \  *  / V  n_d  //  V"  )  (16) 


■bars  ARlo  and  AltL^  are  the  average  run  lengths  in 
the  in-control  and  out-of-control  states,  respec¬ 
tively,  and  E(T)'  is  given  by  equation  (2) .  Recall 
that  a  and  1-8  are  functions  of  the  sample  size  n, 
the  control  limit  factor- k,  and  the  magnitude  of 
the  shift.  Kota  that  this  la  a  wulticr iter ion  opti¬ 
misation  problem,  and  would  generally  be  herder  to 
eolve  then  the  optimization  problems  usually  asso¬ 
ciated  with  Che  economic  design  of  control  charts. 
However,  it  avoids  explicit  estimation  of  the  seerch 
costs  ij  and  >3,  and  the  penalty  cost  for  production 
in  the  out-of-control  state  eg,  which  many  analysts 
(lad  difficult. 

toother  possible  formulation  la  to  minimise  the 
sampling  cost  on  a  per  sample  basis  while  simulta¬ 
neously  mexlmislng  ARLg  and  minimising  ARt,.  Thus, 
<11)  becomes 

minimise  t"  -  Sj  ♦  a2n  (14) 

for  all,  and  (12)  end  (13)  ere  unchanged.  This 
formulation  of  the  problan  involves  only  n  and  k. 


whara  6  la  the  lot  fraction  defective.  The  plot  of 
P(8)  versus  8  for  0  s  8  s  1  is  callad  the  operating 
characteristic  curve  (or  OC  curve)  of  the  plan. 
Traditionally,  we  design  acceptance  sampling  plana 
(l.e.,  choose  n  and  c)  so  that  P(6)  passes  through 
or  near  two  points,  such  as  the  familiar  producer's 
and  consvntcr's  risk  points  (AQL.a)  end  (LTPD.B) . 

When  N  Is  lsrgs  relative  to  n,  the  hypergeometric 
distribution  in  (16)  esn  be  repleced  by  the  blnosulal. 

Prior  Distributions  tor  Attribute  Sampling 

Lot  fraction  defective  8  is  e  function  of  two 
sources  of  variability;  the  variability  of  p  the 
process  fraction  defective,  end  the  variability  of 
8  about  p.  It  is  usually  convenient  to  assume  that 
lot  quality  has  a  mixed  binomlel  distribution;  that 
is,  each  lot  is  produced  by  s  production  process 
that  is  la-control  at  the  level  p,  but  p  varies  from 
lot  to  lot  according  to  a  probability  distribution 
f(p).  the  distribution  f(p)  is  often  called  the 
prior  distribution  for  p  cr  the  process  curve.  It 
Is  sxtrmnely  important  to  note  chat  if  the  process 
quality  is  stable  such  that  f(p)  *  1  whan  p  ■  po 


and  f  <p)  ”  0  when  p  y  p,-v ,  then  * here  is  no  need  for 
sampling  (this  result  is  called  Mood's  theorem;  see 
*>od  [15]). 


where  a  is  a  constant  proportional  to  the  variable 

cost  of  sampling  and 


While  there  are  winy  possible  choices  tor  the 
prior  distribution  for  p,  some  of  the  Bora  important 
are  the  following:  the  continuous  beta  distribution. 


where  p  *  *(-u) .  In  practice,  it  Is  important  to 
know  how  accurately  the  prior  distribution  must  be 
specified.  Generally,  the  analyst  does  not  possess 
sufficient  information  about  the  process  to  specify 
the  nr lor  with  great  confidence.  Fortunately,  most 
results  indicate  that  precise  specification  of  the 
prior  la  not  critical,  provided  that  a  reasonable 
distribution  la  chosen.  Continuous  prior  distribu¬ 
tions  are  generally  thought  to  ba  more  appropriate 
than  discrete  ones,  and  the  beta  distribution  (17) 
has  been  used  extensively.  However,  when  the  under¬ 
lying  quality  characteristic  Is  a  continuous  variable 
that  la  normally  dlstrlbutad  within  each  lot  and  the 
mean'  of  thia  quality  characteristic  also  haa  a  nor¬ 
mal  prior,  then  the  prior  distribution  for  p  has 
the  form  (19) .  The  beta  and  normal-generated  dis¬ 
tributions  can  have  vary  different  shapes  for  the 
earns  mean  and  variance,  and  so  significant  differ¬ 
ences  in  the  optimal  sampling  plana  may  result.  For 
further  discussion  of  prior  distributions,  see 
WatharlU  sad  Chiu  [21]. 


The  Mi  lor  A 
Attributes 


aroachss  to  Simla  Sampling 


Given  a  suitable  prior  distribution  and  a  sat 
of  coats  or  losses  associated  with  sampling  plan 
operation.  It  la  desirable  to  choose  the  sampling 
plan  parMstara  that  minimize  the  totel  cost.  Per¬ 
haps  the  most  widely  used  and  detailed  model  la 
that  of  Guthrie  and  Johns  [9].  A  simplified  form 
of  this  modal  la  also  presented  by  Held  [11].  the 
modal  Is  a  linear  cost  model.  All  linear  coat  for¬ 
mulations  lead  to  tha  sane  aspect ad  loss  function. 
Hald  utilizes  the  concept  of  break-even  Quality  pr, 
a  fraction  daf active  value  at  which  It  la  just  as 
costly  to  secapt  as  to  rsjact  tha  lot.  Tha  expected 
loss  par  lot  la 


* n 


1  “  an  +  (W-n)i  /  (p_-p)[l-FCp)3f(p)dp 


*/ 
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is  the  probability  of  accepting  a  lot  of  quality  p. ’  - 
Miniraiting  1  with  respect  to  n  and  c  will  produce 
the  optimal  sampling  plan.  This  is  sometimes  called 
the  Bayesian  approach  to  designing  a  sampling  plan. 

The  major  focus  of  the  research  In  this  area  has 
started  with  the  Guthrie  and  Johns  model.  Hald  haa 
bean  a  major  contributor  in  the  field,  along  with 
eome  of  his  coworkers.  The  major  emphasis  has  been 
on  finding  asymptotic  relationships  between  n  and  c 
for  various  process  curves,  and  in  producing  tables 
suitable  for  use  by  professional  practitioners.  Hald 
[10]  gives  a  number  of  significant  findings  in  his 
I960  paper.  One  part  of  the  paper  Investigates  the 
compound  hypergeometrlc  distribution;  that  la,  tha 
probability  distribution  of  the  number  of  defectives 
d  in  a  random  sample  of  size  n  given  a  prior  distri¬ 
ct  Ion.  The  second  part  of  the  paper  is  directed 
towards  actually  finding  optimum  sampling  plana  for 
rectangular,  beta,  and  double  binomial  priors.  A 
general  solution  is  given  assuming  the  linear  cost 
model  and  inequalities  are  given  for  n  and  c.  In 
1965,  Hald  [11]  provided  tables  for  the  double  bino¬ 
mial  prior  and  In  1968,  Hald  [12]  provided  tables 
for  tha  beta  prior,  along  with  asymptotic  relation¬ 
ships  between  n  end  c,  and  between  8  and  n. 

Despite  tha  significance  of  Held's  work,  the 
tables  he  haa  provided  are  often  difficult  Co  use 
because  of  the  large  amount  of  Information  required. 
An  alternative  approach  consists  of  formulating  an 
appropriate  coat  model  and  optimizing  it  for  a  spe¬ 
cific  problem  using  direct  search  methods.  Consi¬ 
derable  work  In  this  area  has  been  dona  by  G.K.  Ben¬ 
nett,  K.E.  Casa,  and  J.W.  Schmidt  and  their  students. 
Tor  example,  their  1972  and  1975  papers  [3] [17]  deve¬ 
lop  economic  models  for  single  sampling  plans  for 
dealing  simultaneously  with  multiple  attributes.  Tha 
coat  modal*  consist  of  a  component  representing  In¬ 
spection  costa,  a  component  representing  the  expec¬ 
ted  cost  of  lot  rejection,  and  a  component  represen¬ 
ting  the  Expected  coat  of  lot  acceptance.  Pattern 
search  la  used  for  model  optimization.  In  general, 
direct  search  method*  are  a  very  affective  approach 
for  determining  economically  optimal  acceptance 
sampling  plane. 

The  economic  impact  of  tha  disposition  policy  for 
rejected  lot*  has  alao  been  investigated.  While 
there  are  a  nuabar  of  possible  lot  disposition  pol¬ 
icies,  the  two  cases  that  have  been  Investigated  most 
extensively  are  where. rejected  lots  are  either 
scrapped  or  screened  (100  percent  Inspection).  For 
work  on  lot  disposition  policies,  a**  [18]  and  [20]. 

Other  Work 

Thera  Hess  bean  many  other  studies  devoted  to 
tha  economic  design  of  acceptance  sampling  plane; 
for  Instance,  a  1975  survey  paper  by  Watharill  and 
Chiu  [21]  cltae  25)  references.  While  most  of  tha 
work  focuee*  on  tingle  sampling  for  attribute*,  that* 


has  been  lose  research  on  the  economic  design  of 
■ore  sophisticated  sampling  plans.  Foe  example, 
Stewart,  Montgomery  and  Rallies  [20]  describe  a  proce¬ 
dure  for  the  selection  of  double  sampling  plans  for 
attributes  based  on  prior  distributions  and  costs. 
Models  are  presented  for  the  esses  where  rejected 
lots  are  either  screened  or  scrapped.  They  note  that 
there  is  often  little  difference  In  cost  between 
economically  optimal  single  and  double-sampling  plans. 
However,  when  sampling  ccst3  are  large,  double  sam¬ 
pling  plans  have  much  to  offer.  They  also  observe 
that  arbitrary  double-sampling  plans,  such  as  those 
in  MIL  STD  105D,  may  bo  very  far  from  economically 
optimum. 

The  effect  of  Inspection  error  on  sampling 
plan  design  has  also  received  considerable  attention. 
If  an  Inspector  misclassif lee  good  and  bad  items  with 
constant  probabilities,  the  effect  Is  to  translate 
the  OC  curve  of  the  eampllng  plan,  so  that  the  actual 
or  effective  OC  curve  la  somewhat  different  from  the 
nominal  or  advert lead  OC  curve.  If  the  probabilities 
of  mlsclaealf icatlon  are  known,  then  one  may  directly 
Incorporate  this  Information  into  the  economic  design 
of  the  sampling  plan.  Generally,  the  presence  of 
Inspection  errors  implies  that  larger  samples  are 
necessary.  A  good  review  of  the  literature  In  this 
srsa  Is  in  Dorris  and  Foote  16]. 

Vary  little  attention  has  been  given  to  accep¬ 
tance  sampling  by  variables.  Variables  sampling  is 
.not  as  widely  used  in  practice  as  attributes  sampling 
for  several  reasons: 

1.  Variables  measurement  la  often  more  expensive 
(difficult)  than  attributes  measurement. 

2.  A.  asperate  plan  must  be  used  for  each  quality 
characteristic . 

3.  The  estimation  of  fraction  defective  assumes  a 
normally  distributed  quality  characteristic.  If 

,  this  assumption  is  violated,  the  tell  areas  may 

!  he  dramatically  affected  and  the  resulting  frac¬ 

tion  defective  estimate  grossly  In  error. 

A.  It  la  possible  to  reject  a  lot  without  actually 
finding  any  defectives,  end  this  often  upsets 
both  producers  and  consumers. 

However,  remember  that  variables  sampling  can  greatly 
reduce  the  required  sample  else,  and  that  it  does 
generally  provide  better  information  about  the  lot  . 
or  process  quality.  For  work  on  the  economic  design 
of  variables  eampllng  plan,  see  Allot,  Schmidt  and 
Bennett  [1],  Schmidt,  Bennett  and  Case  [18],  and 
Schmidt,  Case  and  Bennett  [19].  Reference  [1]  deals 
with  the  situation  where  the  quality  character  let ice 
•re  a  mixture  of  attributes  and  variables.  In  ell 
of  these  studies,  the  approach  taken  is  to  formulate 
a  cost  modal  and  optimise  it  via  direct  search 
methods. 


Dm  last  20  years  have  seen  the  development  of 
Burner one  techniques  tor  the  design  of  process  control 
and  acceptance  sampling  schemes  based  on  economic 
considerations.  However,  the  indication  la  that  very 


very  few  pr.-.ct  It  loners  have  Implemented  any  of  these 
techniques  (for  example,  see  the  survey*  by  Sanlgs 
and  Shlrland  [16]  and  Chiu  and  Uetherlll  [5]).  This 
la  surprising,  as  most  quality  assurance  managers 
claim  that  cost  reduction  and  increased  productivity 
is  a  major  objective  of  their  function.  In  many 
cases,  an  experienced  engineer  could  design  an  appro¬ 
priate  technique,  perhaps  even  one  that  is  nearly  - 
economically  optimal,  but  the  use  of  a  formal  econo¬ 
mic  model  to  assist  the  analyst  is  a  much  more  pre¬ 
cise  approach,  leaving  lees  to  judgment.  Further¬ 
more,  there  Is  often  a  significant  economic  panalty 
associated  with  the  “standard",  judgment  designs  so 
frequently  used  In  practice. 

The  Implementation  of  these  techniques  requires 
a  computer  program  of  the  cost  model  and  the  optimi¬ 
zation  procedure.  The  lack  of  availability  of  suit¬ 
able  computer  software  has  certainly  slowed  practical 
iapl ament at Ion.  This  is  an  Important  gap  between 
theory  and  practice  that  must  be  filled  before  the 
economic  design  of  quality  assurance  techniques  will 
baccate  widespread. 

Many  practitioners  sre  reluctant  to  use  these 
techniques  because  of  the  difficulty  In  estimating 
costs,  prior  distributions,  and  other  model  parame¬ 
ters.  Fortunately,  costs  sod  prior  distributions  do 
not  have  to  be  estimated  with  high  precision,  although 
some  other  parameters,  such  as  the  magnitude  of  the 
process  shift,  require  more  careful  determination. 
Sensitivity  analysis  of  the  specific  model  could 
help  the  analyst  discover  which  parameters  ere  cri¬ 
tical  In  his  specific  application.  The  availability 
of  efficient,  interactive  computer  software  would  be 
of  significant  value  In  this  respect. 

.  Me  will  illuetzate  the  use  of  e  simple,  inter¬ 
active  computer  program  for  the  optimal  economic 
design  of  en  x-chert.  The  program  assumes  that  the 
Duncan  single-cause  process  modal  is  appropriate, 
and  requests  the  user  to  Input  values  of  ay,  aj,  *3, 
a],  %4,  >,  6,  g,  and  D.  The  optimal  control  limit 
width  k.  Interval  between  samples,  and  minimum  cost 
per  unit  time,  are  calculated  and  displayed  for  a 
range  of  sample  sizes,  along  with  the  a-rlsk  and 
power  of  tne  control  chart  for  each  (n,k,h)  combina¬ 
tion.  The  economically  cptlmal  control  chart  design 
may  be  found  by  inspection  of  the  coat  function  val¬ 
ues  to  find  the  minimum.  The  output  provided  enables 
the  analyst  to  determine  the  sensitivity  of  the  cost 
surface  in  the  vicinity  of  the  optimum.  The  program 
Is  vrittsn  In  FORTRAN  for  b  CDC  CTBEK-74  computer, 
and  involves  less  than  100  lines  of  code.  Execution 
times  are  typically  one  second  or  lees.  Further 
details  of  the  program  sre  ivailable  from  the  author. 

Consider  s  manufacturer  of  non-returnable 
glass  bottles  for  packaging  a  carbonated  soft  drink 
beverage.  The  wall  thickness  of  the  bottles  tr  an 
Important  quality  characteristic.  If  the  well  is 
too  thin,  internal  pressure  generated  during  filling 
will  causa  the  bottle  to  burst.  The  manufacturer 
has  bean  using  5c  and  R  charts  with  n»5,  k*3,  and 
taking  staples  every  h»2  hours  fo  control  the  pro¬ 
cess.  He  wishes  to  compare  this  with  an  econcmically- 
crpt  1ms  1  design  and  estimate  the  savings. 

Based  an  an  analysis  of  ths  quality  control 


technician*'  salaries,  tho  costs  of  the  test  equlp- 
nent ,  and  the  length  of  time  required,  it  in  esti-a- 
ted  that  the  fixed  cost  of  taking  a  sample  is  Si. 00, 
and  the  variable  cost  of  sampling  is.  approximately 
$0.10  pet  bottle.  It  takes  about  one  minute  (0.0167 
hours)  to  measure  and  record  the  wall  thickness  of  a 
bottle. 

The  process  is  subject  to  several  typos  of 
assignable  causes  resulting  in  large  shifts,  typi¬ 
cally  of  about  two  standard  deviations.  (This  is 
why  the  "standard"  x-chart  design  n»5,  k*3  has  been 
used  -  it  is  known  to  be  reasonably  effective  for 
large  shifts).  Ihe  previous  six  months  of  line  per¬ 
formance  data  Indicate  that  the  mean  tine  between 
process  breakdowns  is  about  20  hours.  Thus  \«0.05  is 
the  parameter  of  the  exponential  distribution.  The 
average  rime  to  investigate  an  action  signal  is  one 
hour.  Real  alarms  Incur  a  $25.00  search  cost,  on 
tho  overage,  while  falsa  alarms,  which  ara  more  dif¬ 
ficult  to  check  out,  incur  a  search  cost  of  $50.00. 

The  soft  drink  bottler  to  whom  the  bottlas  are 
sold  has  a  policy  of  backcharging  the  bottle  manu¬ 
facturer  for  the  costs  of  cleanup  and  lost  production 
when  an  excessive  number  of  defective  bottles  burst 
duping  filling.  The  past  six  months'  experience 
indicates  that  $100  per  hour  is  a  reasonable  estimate 
of  the  penalty  coat  of  operating  in  Che  out-of -control 
•tats. 


Figure  !  snows  the  ertual  computer  input  and  out¬ 
put  for  this  example.  Note  that  the  optlimim  design 
has  ty- 5,  k*2.0v,  h-0.76  hours  (about  45  minutes),  and 
■  -Inimum  cost  of  $10.38  per  hnwr.  Thus  the  implica¬ 
tion  is  that  the  manufw<<-urer's  existing  control  chart 
v’sign  is  good  wl*h  respect  to  n  and  k,  but  that  he 
need*  to  samp.  =  more  frequently.  The  actual  cost  of 
his  currant  procedure  is  $12.05  per  hour,  or  about  'a 
16  percent  penalty  cost.  While  the  savings  per  hour 
seem  small,  this  is  a  continuous  production  process 
operating  three  shifts  per  day,  and  the  annual  sav¬ 
ings  are  over  $14,000.00. 

■  • 

After  looking  at  the  optimal  x-chart  design,  the 
bottle  manufacturer  suspects  that  he  may  have  Incor¬ 
rectly  estimated  the  penalty  cost  of  out-of-control  . 
production  (a <,),  and,  at  worst,  he  may  have  underes¬ 
timated  this  parameter  by  50  percent.  Therefore,  he 
reruns  the  program  with  a&*$150.00  to  investigate 
the  affect  of  mlsepecifylng  this  parameter.  The  com¬ 
puter  output,  shown  in  Figure  2,  indicafes  that  the 
optimum  design  it  now  n>5,  V-2.99,  h-0.62,  with  a 
minimum  cost  of  $13.88  per  hour.  Note  that  the  pri¬ 
mary  effect  of  lncreaaing  aq  by  50  parcent  is  to  de¬ 
crease  the  interval  between  samples  from  45  minutes 
to  37  minutes.  The  program  could  be  used  in  this 
manner  to  quickly  and  easily  investigate  the  effects 
of  errora  in  specifying  any  of  the  model  parameters. 
Based  on  this  analysis,  the  manufacturer  electa  to 
adopt  a  45  minute  sampling  interval,  because  It  is 
administratively  simple  and  it  will  not  require  any 
additional  quality  control  technicians. 


1ENTER  SYSTEM  P ARAMETE  RSSAl,A2r A3  r  A3PR  Z ME  r  A4  r  LAMBDA  »  BELT  A  »  G  r  D 


1. 0*0.1 >23. 0*50. Or  100. Or 0. 05r 2. 0r0.0167r 1.0 


N 

OPTIMUM  K 

OPTIMUM  H 

ALPHA 

POWER 

COST 

1 

2.30 

.43 

.0214 

.3821 

14,71 

2 

2.32 

.57 

.0117 

.6211 

11.91 

3 

2.48 

.66 

.0074 

.7835 

10.90 

'4 

2.84 

.71 

.0045 

.8770 

10.31 

S 

2.99 

.76 

.0028 

.9308 

10.38 

4 

3.13 

.79 

.0017 

.9616 

10.39 

7 

3.27 

.82 

.0011 

.9784 

10.48 

0 

3.40 

.85 

.0007 

•  98S0 

10.60 

t 

3.53 

.87 

.0004 

.9932 

10.75 

10 

3.64 

.89 

.0003 

.9961 

40.90 

11 

3.78 

.92 

.0002 

.9978 

11.06 

12 

3.90 

.94 

.0001 

.9988 

11.23 

13 

4.02 

.96 

.0001 

.9993 

11.39 

14 

4.14 

.98 

.0000 

.9996 

11.56 

IS 

4. 29 

1.00 

.0000 

.9998 

11.72 

.943  CP  SECONDS  EXECUTION  TIRE 


Figure  1.  Sample  Computer  Output 


IENTER  SYSTEM  PARAMETERS: A1 »A2 , A3» A3PRIME ,AA, LAMBDA . DELTA ,0,D 
? 


.0»0 

.1.25.0. 

50,0. 150. 0.0. 

05.2.0. 

0.0167.1 

.0 

N 

OPTIMUM 

K  OPTIMUM  H 

ALPHA 

POUER 

COST 

1 

2.31 

.37 

.0209 

.3783 

19.1? 

2 

2.52 

.46 

.0117 

.6211 

15.71 

3 

2.68 

.54 

.0074 

.7835 

14.48 

4 

2.84 

.58 

.0045 

.8770 

14.01 

5 

2.99 

.62 

.0028 

.9308 

13.88 

6 

3.13 

.65 

.0017 

.9616 

13.91 

7 

3.27 

.67 

.0011 

.9784 

14.04 

B 

3  •  40 

.69 

.0007 

.9880 

14.21 

9 

3,53 

.71 

.0004 

.9932 

14.41 

10 

3.66 

.73 

.0003 

.9961 

14.62 

11 

3.78 

.75 

.0002 

.9978 

14.84 

12 

3.90 

.77 

.0001 

.9988 

IS. 06 

13 

4.02 

.78 

.0001 

.9993 

15.28 

14 

4.14 

.80 

.0000 

.9996 

15.50 

.417  CP  SECONDS  EXECUTION  TIME 


Figure  2.  Sample  Computer  Output 
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